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ABSTRACT 

a   «n  Tt?  obJectives  of  this  thesis  are  to  1)  determine  the 
deflections  of  the  horizontal  and  vertical  tank  supports 
when  the  ship  is  subjected  to  various  loading  eond???™* 

oJatformf  fo^  **?  3^S5eS    ln  ^  ^™*wZ%^' 
SiSeeiSlv  2h.  S!??  ^ding  $onditio*s.   This  information, 
especially  the  deflections  of  the  support  platform   iq 
needed  to  determine  the  forces  that  Sill  be  tSSSi+JL  * 
the  ship  through  the  articulated  vesical  suppo^svstem  IT 
the  large  spherical  pressure  vessels  that  ca??y  thfcarL 
XSrSS?"  WU1  P6rmit  a  —  o-pleLrarLiyhesinf  ihe 


TrRq  S;^faPer  details  how  the  beam  element  capability  of 
ICES  STRUDL  was  used  to  determine  the  stresses  and  ^flec- 
tions. A  quarter  tank  section  of  the  ship  is  modeled  as  a 
three  dimensional  arrangement  of  beams.  This  model  is  then 
subjected  to  loads  that  simulate  the  loads  the  actual  m«£ 
ter  tank  section  would  experience.  qUar" 


ti™c  IS  ^sults  mficate  that  the  maximum  vertical  deflec- 
?i2S  +S!     support  platform  out  of  a  horizontal  plane  are 
less  than  a  centimeter  for  the  loadings  investigated   The 
circular  hold  deforms  out  of  its  circular  shSpe^n  ail  thS 

ment  of  nil  °™   W°+ld  ?Xpect*   The  maximum  horizontal  move- 
ment of  any  support  point  on  the  circular  suDDort  nia+f^! 

slresleVa^  JaSX^S™-  .R%lativ"y"i^PL°n|i?idaInaim 
cent??nne  of  ?h.  +^\theATnain  deck  at  the  transverse 
£«?  +  i ine  °f  *he  tank-   Areas  of  high  bending  stresses  are 
pointed  out  and  should  be  further  investigated  usinsmeans 
other  than  a  beam  element  model  due  to  modeling  difficSuLs. 
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INTRODUCTION 


The  first  carriage  of  liquified  natural  gas  (LNG)  oc- 
curred in  1958  on  an  experimental  basis  with  the  conversion 
of  an  old  CI  dry  cargo  ship  into  the  Methane  Pioneer.   Suc- 
cessful commercial  ventures  took  place  in  1964  when  the 
first  new  ships  specifically  designed  for  the  transporta- 
tion of  LNG  were  completed  and  put  into  service.  At  the 
present  time  there  are  approximately  20  LNG  tankers  in  use. 

The  reason  for  this  large  growth  of  interest  in  the 
marine  transportation  of  liquified  natural  gas  has  been  the 
world  wide  shortages  that  have  been  caused  not  only  by  in- 
creasing general  demand,  but  also  because  of  the  demand  for 
"clean"  fuels  which  is  inspired  by  the  ecological  concern 
being  expressed  at  the  present  time.   Natural  gas  is  the 
cleanest  burning  fuel  that  is  available  at  the  present  time 
with  the  exception  of  hydrogen. 

Natural  gas  is  often  a  by-product  of  oil  drilling 
operations  and  has  often  been  considered  a  waste  product  in 
the  past,  especially  in  areas  remote  from  consumers.  Trans- 
portation of  this  fuel  in  the  gaseous  state  would  be  eco- 
nomically infeasible  due  to  the  large  volume  required.  How- 
ever, it  was  determined  that  if  the  gas  could  be  liquified, 
the  required  volume  per  unit  weight  could  be  greatly  re- 
duced. 

The  principal  constituent  of  natural  gas  is  methane 
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which  has  a  boiling  point  of  -258.6  degrees  Fahrenheit  and 
a  gas  to  liquid  volume  ratio  of  6301 1  (Ref.  1).   Thus,  if 
the  temperature  of  the  natural  gas  is  reduced  below  -260 
degrees  Fahrenheit,  the  same  amount  of  fuel  could  be  carried 
in  I/630  of  the  volume  of  the  gas  at  room  temperature. 

These  low  temperatures  caused  several  technical  prob- 
lems of  which  the  principal  one  was  probably  brittle  frac- 
ture in  steel.   Thus,  the  tank  material  must  have  a  ductile 
to  brittle  transition  temperature  well  below  -260  degrees 
Fahrenheit.   The  materials  used  for  the  construction  of  the 
LNG  tanks  to  date  have  been  aluminum,  stainless  steel,  Invar 
(36  per  cent  nickel-iron  alloy)  and  9  per  cent  nickel  steel 
(Ref.  12).   The  problem  of  brittle  fracture  required  these 
tanks  to  be  thermally  insulated  from  the  hull  of  the  ship 
which  is  constructed  of  mild  steel. 

There  have  evolved  two  basic  tank  types  by  which  LNG 
is  carried — the  self-supporting  tank  and  the  membrane  tank. 
There  are  several  differing  designs  of  each  basic  tank  with 
varying  support  systems,  insulations,  tank  shapes,  materi- 
als, and  type  of  secondary  barrier,  if  any. 

Membrane  tanks  are  designed  such  that  they  are  liquid 
tight  only,  with  the  intention  that  the  normal  loading 
forces  would  be  transmitted  through  the  thermal  insulation 
to  the  hull  structure  on  the  other  side.  Membrane  tanks 
require  a  secondary  barrier.   The  secondary  barrier  must  be 
designed  to  contain  a  catastrophic  failure  of  the  membrane 
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for  a  reasonable  amount  of  time  (Ref.  <?)• 

There  are  two  basic  types  of  self-supporting  tanks — 
pressure  vessels  and  free  standing  tanks.   Pressure  vessels, 
due  to  their  shapes,  are  amenable  to  analytical  stress  a- 
nalysis  and  are  constructed  in  accordance  with  pressure 
vessel  requirements.   Free  standing  tanks  are  basically 
tanks,  usually  of  rectangular  or  trapezoidal  shape,  that  are 
set  inside  the  ship's  hull  on  supports.   The  advantage  of 
the  free  standing  tank  over  the  pressure  vessel  is  that  it 
can  conform  more  nearly  to  the  ship's  hull  configuration. 
The  disadvantage  is  that  the  free  standing  tanks  require  a 
secondary  barrier  and  the  pressure  vessels  do  not.  As  a  re- 
sult of  economic  considerations,  some  recent  designs  have 
been  of  very  large  pressure  vessels  having  low  vapor  pres- 
sures with  the  differences  between  the  two  types  of  free 
standing  tanks  becoming  indistinguishable.  The  distinction 
between  free  standing  tanks  and  pressure  vessels  varies 
among  the  regulatory  agencies,  however,  the  use  of  a  vapor 
pressure  of  10  psig  as  a  demarcation  limit  is  common. 

The  United  States  Coast  Guard,  which  is  charged  with 
the  regulation  of  the  marine  industry  from  the  safety  as- 
pect, has  recently  defined  four  types  of  self-supporting 
tanks  (Ref.  16). 

1 .   Independent  Pressure  Vessel  Tanks  —  These  are 

pressure  vessels  with  the  vapor  pressure  P  great- 
er than  EH-20  where  D  is  the  diameter  in  feet  and 
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where  P  is  at  least  40  psig. 

2.  Semi -Independent  Pressure  Vessels  —  This  category 
includes  pressure  vessel  tanks  where  the  vapor 
pressure  is  less  than  for  the  independent  pressure 
vessels. 

3.  Independent  Structurally  Indeterminate  Self -Sup- 
porting Gravity  Tanks  —  These  are  free  standing 

tanks  in  which  the  maximum  vapor  pressure  is  less 

i 

than  10  psig. 

^.   Independent  Structurally  Semi-Determinate  Self- 
Supporting  Gravity  Tanks  —  These  are  large  tanks 
with  vapor  pressures  under  10  psig  that  are  usually 
designed  using  finite  element  and/or  fracture 
mechanic  techniques  to  pressure  vessel  standards. 
This  class  of  tank  generally  will  be  required  to 
have  only  a  partial  liquid  tight,  splash  tight 
secondary  barrier. 

The  proposed  tanks  for  the  Technigaz  design  that  is 
being  analyzed  fall  into  category  4. 
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DESCRIPTION  OF  THE  SHIP  AND  CARGO  TANKS 
» 
The  Ship 

This  analysis  was  undertaken  as  part  of  the  design 
process  for  the  Technigaz  125»000  cubic  meter  methane  car- 
rier design.   The  principle  characteristics  of  the  design 
are  as  follows t 

Length  overall  290.6  meters 

Length  between  perpendiculars  275. 0  meters 

Breadth,  molded  ^4.0  meters 

Depth(at  main  deck)  26.0  meters 

Draft  10.97  meters 

Block  coefficient  0.775 

The  cargo  tanks  consist  of  five  spherical  tanks  ver- 
tically supported  at  their  equator  (Fig.  1  and  2),   The 
spaces  for  these  spherical  tanks  occupy  a  large  portion 
of  the  deck  area  necessitating  the  requirment  of  very  heavy 
plating  (^5  millimeters  thick)  at  the  deck  and  sheer  strake 
junction  to  provide  longitudinal  material  in  tension,  stiff- 
ness against  lateral  loads  from  the  sea,  and  resistance  to 
buckling  in  the  sagging  condition. 

Longitudinal  bulkheads,  located  8.9  meters  inboard 
from  the  outer  hull  except  in  way  of  tank  holds,  provide  the 
lateral  subdivision  for  wing  tanks  which  are  used  for  bal- 
last (Fig.  3).   The  double  bottom  is  composed  of  longitu- 
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FIGURE  2 
MIDSHIP  SECTION  -  FRAME  228 
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FIGURE  3 
TYPICAL  DECK  -  QUARTER  TANK  SECTION 


-16- 


dinal  plates  spaced  3 .28  meters  apart.   In  addition  the 
inner  and  outer  bottoms  are  strengthened  by  longitudinal 
stiff eners  spaced  .82  meters  apart. 

Transversely  the  hull  is  stiffened  by  deep  web  frames 
in  the  double  bottom  and  sides.   The  frame  spacing  is  .85 
meters  with  a  web  frame  every  2.55  meters.   There  are  six 
decks  in  addition  to  the  double  bottom  and  the  main  deck. 

The  Cargo  Tanks 

The  cargo  tanks  consist  of  five  spherical  tanks  of 
semi-independent  pressure  vessel  type.  Tank  1  has  an  in- 
side diameter  of  31*6  meters  and  a  capacity  of  16,520  cubic 
meters.   Tanks  2,  3,  4,  and  5  have  an  inside  diameter  of 
37.5  meters  and  a  volume  of  27,610  cubic  meters.  This  re- 
sults in  a  maximum  capacity  of  126,965  cubic  meters.  The 
useful  capacity,  using  a  maximum  filling  ratio  of  98$,  is 
12^,450  cubic  meters. 

The  spherical  tanks  are  supported,  in  the  vertical 
direction,  at  the  equator  by  a  patented  system  of  articu- 
lated parallelograms  of  rods  and  arms  (Fig.  4).   This  sys- 
tem is  attached  to  brackets  at  the  tank  equator  and  then 
extend  down  to  the  support  platform.   Horizontal  restraint 
to  pitch  and  roll  is  provided  by  keys  and  keyways.  These 
keys  are  installed  on  two  circles  parallel  to  the  equator 
at  the  approximate  positions  of  the  Tropic  of  Cancer  and 
Tropic  of  Capricorn. 
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FIGURE  4 
TECHNIGAZ  ARTICULATED  SUPPORT  SYSTEM 
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The  purpose  of  the  vertical  support  system  of  articu- 
lated parallelograms  is  to  allow  the  cargo  and  tank  weight 
to  be  supported  without  a  bending  moment  being  introduced 
in  the  shell.   It  also  allows  for  the  thermal  contraction 
of  the  sphere  while  carrying  the  cold  cargo.   The  keys 
change  the  horizontal  forces  due  to  pitch  and  roll  to  tan- 
gential forces  on  the  tank  shell.   The  purpose  of  these  e- 
laborate  support  systems  is  to  allow  only  forces  that  can 
be  structurally  analyzed  to  be  imparted  to  the  spherical 
tanks.   This  allows  the  tank  plating  thickness  to  be  re- 
duced to  the  minimum  possible  thickness  and  also  allows  a 
possible  reduction  in  the  secondary  barrier. 

The  present  proposed  design  scheme  is  to  fit  insulation 
around  the  tank  inside  the  inner  hull  with  the  inner  face  of 
the  insulation  acting  as  a  spray  shield  in  case  of  tank 
leakage.  A  drip  pan  is  fitted  in  the  bottom  of  the  hold  to 
collect  any  drippings  and  to  act  as  a  local  secondary  bar- 
rier. 

One  of  the  purposes  of  this  thesis  is  to  determine  the 
deflections  of  the  vertical  support  platform  and  the  key- 
ways  in  order  that  the  forces  imparted  to  the  sphere  may  be 
determined  prior  to  a  structural  analysis  of  the  shpere, 
which  is  required  before  a  reduction  in  the  secondary  bar- 
rier will  be  permitted  by  the  regulatory  agencies. 
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MATHEMATICAL  MODEL 

General  Discussion 

The  solution  method  that  appeared  particularly  suited 
to  the  task  of  solving  for  the  deflections  of  the  horizon- 
tal and  vertical  tank  supports  as  well  as  solving  for  the 
stress  levels  in  the  various  structural  members  was  some 
type  of  beam  or  plate  element  solution.   Recent  advances 
and  knowledge  in  this  field  have  been  extensive.   Predicted 
stress  levels  obtained  by  this  type  of  analysis  have  shown 
a  good  correlation  with  actual  results  obtained  from  strain 
gages  on  ships  at  sea. 

It  is  fortunate  that  these  elemental  solutions  have 
been  developed  since  the  traditional  naval  architectural 
methods  of  stress  calculations  would  probably  have  been 
inadequate  for  the  unconventional  geometry  of  this  vessel. 
The  fact  that  the  section  modulus  varies  with  length  (large 
portions  of  the  decks  are  cut  out  to  make  way  for  the  spher- 
ical tanks)  makes  standard  calculations  difficult  if  not 
impossible. 

There  are  several  element  type  computer  programs  avail- 
able that  have  either  been  developed  specifically  for  ships 
or  that  are  general  purpose.   Some  of  these,  such  as  ICES 
STRUDL-II,  developed  by  Massachusetts  Institute  of  Technol- 
ogy for  use  on  IBM  360/370  (Ref.  7  and  8)  and  DAISY  (Ref. 
10),  developed  by  the  University  of  Arizona  in  conjunction 
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with  American  Bureau  of  Shipping  for  use  on  both  CDC  6000 
series  and  for  the  UNIVAC  1108  computers,  have  user  oriented 
languages.   Other  programs  exist,  especially  in  the  aero- 
space field,  where  the  stiffness  matrix  is  assembled  by 
hand  and  a  computer  is  used  for  the  solution. 

The  decision  was  made  to  use  the  M.I.T.  developed  ICES 
STRUDL-II  (Integrated  Civil  Engineering  System-Structural 
Design  Language)  program  which  is  a  series  of  computer  pro- 
grams that  can  be  used  in  conjunction  with  each  other  to 
solve  problems  in  structural  engineering.  The  STRUDL  pro- 
gram was  selected  primarily  because  of  time  consideration. 
The  time  required  to  develop  a  finite  element  program  par- 
ticularly suited  to  this  ship  structure,  would  have  been 
prohibitive.  Another  advantage  of  STRUDL  is  that  the  lan- 
guage is  easily  understandable  to  an  engineer  in  that  the 
words  and  phrases  in  the  input  are  such  that  their  meaning 
is  usually  self-evident. 

Several  different  types  of  analysis  procedures  are 
currently  available  in  STRUDL  for  solving  framed  structures 
and  continuum  mechanics  problems.   The  member  stiffness 
matrix  in  the  frame  analysis  is  computed  from  beam  theory, 
while  continuous  mechanic  problems  are  solved  with  the  fi- 
nite element  capability  and  the  element  stiffness  matrix 
is  computed  from  energy  considerations.   STRUDL  allows  for 
the  mixing  of  members  and  elements  provided  they  have  the 
same  number  of  degrees  of  freedom  per  joint  and  these  free- 
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doms  correspond  in  type  and  direction. 

Analyses  of  ship  structures  have  been  carried  out  using 
beam  elements  as  well  as  finite  elements.   Stiansen  and 
Elbaloute  (Ref,  15)  did  a  thorough  study  of  SL-7  container 
ship  design  using  DAISY.   Fenton  (Ref.  4)  used  the  finite 
element  capability  of  STRUDL  to  obtain  stresses  in  a  cata- 
maran cross  structure.  Zoller  (Ref.  17)  combined  beams  and 
finite  elements  in  a  tanker  study  with  satisfactory  results. 

It  is  very  important  to  use  the  analysis  procedure  that 
is  best  suited  for  the  problem  since  the  number  of  elements 
or  members  required  for  even  a  coarse  node  network  of  a 
complex  structure  can  become  extremely  large  with  the  stiff- 
ness matrix  computer  core  requirements  exceeding  the  capac- 
ity of  the  computer. 

The  decision  was  made  to  use  beam  elements  in  this 
study.   The  main  reason  being  that  the  complexity  of  the 
structure,  caused  by  the  circular  hold  for  the  spherical 
tanks,  required  a  fairly  fine  mesh  with  all  joints  having 
six  degrees  of  freedom  to  get  any  meaningful  deflections. 
This  would  have  required  an  extremely  large  number  of  ele- 
ments if  the  finite  element  capability  was  utilized,  since 
plate  bending  elements  would  have  to  be  superimposed  on 
plane  stress  elements  to  account  for  loads  in  all  three  co- 
ordinates. Any  combination  space  frame  and  finite  element 
would  be  completely  out  of  the  question  for  the  same  reason. 

The  stiffness  analysis  carried  out  by  STRUDL  is  a 


-22- 


linear,  elastic,  static,  small  displacement  analysis.  The 
procedure  used  requires  the  specification  of  the  model  ge- 
ometry s  the  member  and  joint  loads,  the  member  properties 
and  various  applicable  constants,  such  as,  Young 8s  modu- 
lus and  Poisson's  ratio.   The  analysis  used  treats  the  joint 
displacements  as  unknowns.   When  a  space'  frame  is  specified 
as  the  structural  type,  as  it  was  in  this  problem,  there 
are  six  unknowns  at  each  joint — three  displacements  and 
three  rotations. 

The  stiffness  analysis  is  comprised  of  the  following 
steps: 

1 .  Consistency  checks  of  members 

2.  Generation  of  the  stiffness  matrix 

3.  Processing  of  member  loads 

4.  Stiffness  matrix  assembling 

5.  Processing  of  the  joints 

6.  Solution  of  the  matrix 

7.  Joint  displacement  processing 

8.  Processing  of  member  stresses 


-23- 


Beam  Model 

The  area  of  the  midship  section  around  Tank  3  was  the 
section  selected  for  study.   This  area  was  chosen  because 
the  maximum  bending  moment  occurred  in  this  region  for  the 
loading  conditions  that  were  to  be  analyzed.   This  area 
would  then  give  the  maximum  amount  of  deflection  of  the 
support  platform  as  well  as  the  largest  stresses — the  two 
parameters  that  were  to  be  determined. 

In  order  to  keep  the  bandwidth  of  the  stiffness  matrix 
as  small  as  possible  and  still  get  valid  results  it  is  ad- 
vantageous to  limit  the  size  of  the  model.   First  it  was 
possible  to  limit  the  model  size  by  the  use  of  symmetry  of 
both  the  structure  and  the  loading  of  the  ship  about  the 
centerline  of  the  vessel  for  the  cases  being  investigated. 
This  cut  the  size  of  the  model  in  half.   Furthermore,  with 
approximate  symmetry  of  the  moment  curve  and  internal  struc- 
ture of  the  ship  about  frame  228,  it  was  decided  that  suf- 
ficient symmetry  did  exist  to  further  reduce  the  model. 
Thus  the  model  could  be  reduced  to  a  quarter  tank  section 
extending  from  frame  195  to  frame  228  with  the  necessary 
boundary  conditions  at  frame  195  to  simulate  the  remainder 
of  the  vessel. 

All  joints  in  the  transverse  plane  of  symmetry  through 
the  center  of  the  tank  and  in  the  longitudinal  plane  of 
symmetry  along  the  centerline  of  the  vessel  must  be  sup- 
ported.  STRUDL  assumes  that  all  joints  specified  as  sup- 
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port  joints  to  be  rigidly  supported  according  to  structural 
type  with  no  displacements  or  rotation  unless  specifically 
released  in  a  certain  direction.   Thus,  to  insure  an  accu- 
rate representation  of  the  actual  ship,  it  was  necessary  to 
release  the  constraints  in  the  Y  and  Z  direction  as  well  as 
the  rotation  constraint  about  the  X  axis  for  joints  in  the 
transverse  plane  of  symmetry.   In  the  longitudinal  plane  of 
symmetry,  it  was  necessary  to  release  the  force  constraints 
in  the  X  and  Z  direction  in  addition  to  the  rotational  con- 
straint about  the  Y  axis.   The  origin,  the  center  of  the 
tank  hold  in  the  double  bottom,  remained  completely  fixed 
to  provide  a  reference  point. 

The  structural  members  of  the  ship  were  idealized,  as 
indicated  previously,  by  beams.   The  object  of  this  modeling 
process  was  to  accurately  represent  the  ship,  while  at  the 
same  time  trying  to  keep  the  number  of  beam  elements  to  a 
minimum  because  of  computer  space  considerations.  The  model 
that  was  finally  settled  upon  idealized  two  actual  ship 
decks  as  one  deck  in  the  model.   In  the  horizontal  plane 
the  flanges  of  the  beam  element  were  the  deck  plating  and 
the  vertical  members  formed  the  web  of  the  beams.   The 
effective  breadth  of  plating  was  calculated  for  bending  in 
both  the  transverse  and  longitudinal  directions.   Examples 
of  beam  elements  formed  from  the  double  bottom  and  support 
platform  are  shown  in  figures  5  and  6. 

The  vertical  beam  elements  were  handled  in  basically 
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LONGITUDINAL 


FIGURE  5 
DERIVATION  OF  BEAM  ELEMENT  CROSS-SECTIONS 
FROM  CELLULAR  DOUBLE  BOTTOM 
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FIGURE  6 
DERIVATION  OF  BEAM  ELEMENT  CROSS  SECTIONS 
FROM  SUPPORT  PLATFORM  DECKS 
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the  same  way.   Generally  the  "beam  elements  consisted  of  tee 
bars  welded  to  plating.   This  was  true  for  the  outer  verti- 
cal hull,  the  bulkhead  at  the  inner  web  frames,  and  the  cir- 
cular support  platform.   Figure  7  depicts  the  geometric 
shape  of  the  beam  model.   Figures  8  through  16  show  the  in- 
dividual joints  and  member  numbers  for  the  various  joints 
and  bulkheads  of  the  model. 

The  major  problem  that  developed  in  the  modeling  pro- 
cess involved  the  circular  support  platform.   One  charac- 
teristic of  beam  elements  is  that  the  ends  must  always  be 
at  a  joint.   Thus  it  was  necessary  that  the  vertical  mem- 
bers at  least  butt  on  a  longitudinal  or  transverse  member 
of  the  double  bottom.   These  joints  have  numbers  75»  H5i 
156,  198,  251,  and  253  in  the  model(Fig.  8).   This  restric- 
tion prohibited  the  vertical  members  from  being  equally 
spaced  as  in  the  actual  ship.   The  equal  spacing  in  the 
ship  design  was  accomplished  by  the  use  of  vertical  stiff - 
ner  plates  in  the  double  bottom  under  the  vertical  members. 

The  model  did  not  take  into  account  such  details  as 
the  rounded  bilge  or  rounded  deck  edge,  however,  the  effect 
of  these  idealizations  were  judged  to  be  negligible  when 
compared  to  mesh  used. 

It  should  be  pointed  out  here  that  the  point  of  con- 
nection of  the  beam  elements  is  at  the  neutral  axis  of  the 
element.   This  is  the  reason  for  the  small  difference  in 
the  actual  beam  and  depth  when  compared  to  the  model  scale 
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Straight  line  indicates  beam  element 


FIGURE  7 
THREE  DIMENSIONAL  VIEW  OF  THE  BEAM  MODEL 
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FIGURE  8 
BOTTOM  DECK  OF  BEAM  MODEL 
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FIGURE  9 
SECOND  DECK  OF  BEAM  MODEL 
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FIGURE  10 
THIRD  DECK  OF  BEAM  MODEL 
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FIGURE  11 
FOURTH  DECK  OF  BEAM  MODEL 


-33- 


Sol 

yoj 

r°3 

Soy 

ST3JT 

i"o  « 

CO? 

5»8 

© 

<&> 

!E> 

5sj) 

w> 

<&> 

O) 

(lg)         @ 

«3           SIS 

rn 

«r 

£)< 

fl» 

S'8 

5J» 

J"V» 

<K> 

<5e>    na 

$52>      Sa 

(?v§)     S/y 

(fe&  6-/^- 

&*>     «< 

<|$S>     SV7 

<5^>       frh 

56/           1 

sw 

yv* 

**«- 

**7 

SI* 

J9S) 

tf?S) 

g^> 

(S2>         <£h 

r«, 

SJo 

5a  I 

JJZ. 

4s^ 

5*1 

Em* 

5n 

S3 

jtrNy?) 

€?£> 

&h           <&) 

\        »3 

Sit 
ST<t 

51 1 

SS~ « 

N©5) 

>ft»M 

S3)        <SS| 

\    s*' 

5a  8 

«r\ 

sn 

£&       (Sb 

\      r" 

fc*LE    (meTsts) 

5%  \           a> 

0 

i     x     3      *     r 

ys&   fiS) 

1               4 

J 

\         So 

S1*,  1              Slo 

'  * 

1 

_   _P3>  <£? 

»»•'» 


le 


Ml- 


S3, 


4 


(Q)  joint  number 
10  member  number 


FIGURE  12 
MAIN  DECK  OF  BEAM  MODEL 
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FIGURE  13 
OUTER  SKIN  OF  BEAM  MODEL 
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FIGURE  14 
INNER  WEB  BULKHEAD  OF  BEAM  MODEL 
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FIGURE  15 
TRANSVERSE  SECTION  AT  FRAME  195 
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FIGURE  16 
PROJECTION  OF  TANK  HOLD 
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beam  and  depth.   This  point  is  illustrated  by  the  following 
example.   The  depth  of  the  ship  is  26.0  meters.  The  scale 
depth  of  the  model  is  23.95  meters.   This  is  caused  by  the 
fact  that  the  neutral  axis  of  the  double  bottom  beams  are 
1.15  meters  from  the  bottom  plating  and  the  neutral  axis  of 
the  beams  modeling  the  main  deck  are  .90  meters  down  from 
the  main  deck. 

The  location  of  the  average  neutral  axis  was  used  as 
the  location  of  all  the  neutral  axes  of  the  beam  elements 
that  would  be  located  in  a  straight  line  on  the  actual  ship. 
This  can  be  best  illustrated  using  the  outer  web  frames 
where  the  thickness  of  the  shell  plating  changes  as  the 
height  is  increased.   This  change  in  plate  thickness  causes 
the  neutral  axis  of  each  beam  element  to  shift  slightly. 
The  member  properties  of  the  individual  members  were  cal- 
culated using  their  own  neutral  axis.   However,  in  assem- 
bling the  model  geometry,  the  average  neutral  axis  was  used 
allowing  the  members  to  lie  in  a  vertical  straight  line. 

There  are  disadvantages  to  using  the  beam  model  that 
was  outlined  in  the  predeeding  paragraphs.   One  disadvantage 
is  that  the  beam  elements  are  connected  at  the  joints  by  a 
point.   The  actual  area  of  connection  however  may  extend  a 
considerable  distance  from  the  point  that  results  from  the 
intersection  of  the  neutral  axis.   This  is  especially  true 
when  plates  are  being  modeled  as  beams  which  was  the  case 
with  this  model.   The  other  principal  disadvantage  of  this 
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model  was  the  result  of  modeling  two  actual  decks  as  a 
single  deck.   This  would  cause  a  vertical  beam  that  had  a 
support  at  each  end  to  be  modeling  a  member  that  actually 
was  supported  at  three  places  along  its  length.  Along  the 
same  line  a  vertical  member  that  is  subjected  to  a  uniform 
or  linear  force  would  only  be  supported  at  the  ends  while 
the  actual  member  would  be  supported  in  at  least  three 
places.   The  result  of  these  two  disadvantages  would  be  that 
the  local  bending  stresses  in  the  individual  members  would 
be  inaccurate.   Figtsres  17  and  18  illustrate  these  two  dis- 
advantages of  the  beam  model. 

However,  it  should  be  pointed  out  that  these  disadvan- 
tages would  only  affect  local  bending  stresses  to  any  great 
extent.   Overall  results,  such  as  joint  displacements  and 
longitudinal  stresses,  would  be  virtually  unaffected.   In 
addition  the  bending  stresses  that  resulted  from  this  model 
would  give  an  indication  of  areas  that  should  be  looked  at 
more  closely — possibly  with  a  fine  mesh  finite  element 
model  of  that  small  area  using  the  results  of  this  beam 
model  to  determine  the  boundary  conditions. 

The  geometric  model  was  entered  in  the  STRUDL  program 
by  first  numbering  each  joint  and  giving  its  coordinates  in 
the  global  X  Y  Z  space.   The  orientation  of  the  members  was 
determined  by  the  coordinates  of  the  joints  that  were  its 
end  points.   Each  member  has  a  local  coordinate  system 
associated  with  it  for  the  easy  identification  of  input  and 
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junction OF 

TWO  PLATES 


BEAM  MODEL 

CF  ABOVE  FIGURE 


The  light  lines  in  the  top  figure  indicate  material  from 
which  the  member  properties  of  the  beams  were  calculated. 
The  bending  of  the  vertical  plate  about  the  Y  axis  is  not 
accurately  modeled  by  the  two  vertical  beams  since  the 
junction  is  actually  a  straight  line  and  in  the  model  it 
consists  of  two  points. 

FIGURE  17 
MODELING  ILLUSTRATION 
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Figure  (A)  sh 
uniform  load 
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line  showing 
stresses,  the 
not  represent 
of  the  joints 
a  single  beam 
to  two  of  the 


ows  a  beam  supported  at  four  points  with  a 
applied.   Figure  (B)  illustrates  how  the  same 
Id  be  modeled.   As  can  be  seen  by  the  dashed 
the  relative  magnitude  and  sign  of  the  bending 
model  gives  a  local  bending  stress  that  does 
the  actual  case.   However  when  the  deflections 
are  considered  the  modeling  is  accurate  since 
of  the  model  would  have  a  strength  equivalent 
actual  beams. 


FIGURE  18 
MODELING  ILLUSTRATION 
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output.   The  local  positive  x  axis  goes  from  the  starting 
joint  to  the  end  joint.   The  one  remaining  degree  of  free- 
dom for  the  member  was  rotation  about  the  principal  axes. 
This  last  degree  of  freedom  was  specified  by  the  angle  beta. 
The  reference  position  for  the  beta  angle  (beta  =  0)  is 
with  the  local  axis  z  in  the  same  direction  as  the  global  Z 
axis  and  with  the  local  axis  x  parallel  to  the  global  Y 
axis.   The  method  of  determining  the  beta  angle  for  other 
cases  is  outlined  in  Reference  7. 

There  were  many  locations  on  the  actual  ship  where 
additional  stiffners  were  added  to  the  shell  plating  between 
structural  members.   To  account  for  their  effect  the  pris- 
matic cross-sectional  area  of  these  stiffners  were  "blended" 
into  the  plate  to  produce  an  effective  plate  thickness 
(Sample  Calculation  1). 

It  was  not  necessary  to  determine  the  effective  width 
of  plating  for  any  of  the  members  since  no  buckling  of  any 
member  was  anticipated.   However,  it  was  necessary  to  de- 
termine the  effective  breadth  of  the  various  beams,  both  in 
the  y  and  z  directions  in  order  that  a  realistic  moment  of 
inertia  and  section  modulus  could  be  calculated  for  each  of 
the  local  member  coordinates. 

Originally  the  plan  was  to  determine  the  effective 
breadths  using  Schade's  curves  (Ref.  13  and  14).   However, 
the  use  of  these  curves  required  knowledge  of  the  loading 
conditions  prior  to  the  calculation  of  member  properties. 
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The  intention  of  this  work  was  to  do  several  loading  con- 
ditions, (loaded,  ballasted,  hogging,  sagging,  etc.)  and 
it  would  have  "been  completely  infeasible,  time  wise,  to  go 
hack  and  recalculate  the  member  properties  of  all  the  beams 
for  each  loading  condition.   Instead,  the  standard  design 
criteria  for  mild  steel  was  used  where  the  effective  breadth 
of  plating  was  equal  to  sixty  times  the  plate  thickness. 
This  criterion  results  from  the  fact  that  the  effective 
breadth  of  plating  acting  in  association  with  a  structural 
member,  is  influenced  by  the  yield  strength  and  the  modulus 
of  elasticity  of  the  material  (Ref.  11).   This  influence  can 
be  calculated  by  the  formula 

effective  breadth  =  2^JE/F     t 
where  E  =  modulus  of  elasticity  (psi) 
F  =  tensile  yield  strength 

The  above  design  criteria  results  from  the  fact  that  Z-Je/f~ 
equals  60  for  mild  steel. 

STRUDL  requires  that  the  cross-sectional  axis,  the 
torsional  rigidity  coefficient,  and  the  moments  of  inertia 
about  the  local  y  and  z  axes  be  specified  for  the  stiffness 
analysis  of  a  space  frame.   Since  the  sectional  stresses  in 
the  members  were  desired  it  was  also  necessary  to  specify 
the  section  moduli  about  the  local  y  and  z  axes.   The  shear 
areas  in  the  Y  and  Z  directions  were  also  supplied  in  order 
that  the  analysis  would  also  include  shear  deformations. 
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A  computer  program  was  written  to  calculate  the  above 
member  properties  to  avoid  tedious  and  repetitive  hand  cal- 
culations .   The  program  was  capable  of  calculating  the  mem- 
ber section  properties  of  prismatic  beams  of  five  or  fewer 
elements.   The  procedures  used  by  the  computer  program  are 
illustrated  by  an  example  in  Sample  Calculation  2. 

The  beam  element  model  was  constructed  in  such  a  way 
that  it  could  be  expanded  to  handle  non-symmetrical  loadings 
about  the  center,  or  in  other  words,  the  expanded  model 
would  be  capable  of  handling  the  case  when  the  ship  is  in  a 
heeled  position,,   The  final  quarter  tank  section  model  was  a 
space  frame  consisting  of  4l4  beam  members  and  208  joints. 

Information  that  was  requested  as  output  from  the 
STRUDL  program  wasi  (1)  the  forces  acting  on  both  ends  of 
the  member,  (2)  the  reactions  at  support  joints,  and  (3)  the 
displacements  of  each  joint.   In  addition  section  stresses 
of  the  members  were  requested  for  each  end  and  the  center  of 
the  member.   The  section  stresses  were  broken  down  into 
axial  and  bending  stresses  from  which  the  maximum  and  mimi- 
mum  normal  stresses  were  calculated.   The  shear  stresses  at 
the  three  sections  were  calculated  by  dividing  the  shear 
forces  by  the  shear  area. 
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THE  ANALYSIS— HOGGING  CONDITION 

General  Discussion 

The  first  case  to  be  investigated  was  with  the  ship  in 
an  upright  position,  fully  loaded  with  the  wave  crest  amid- 
ships, the  trochoidal  wave  used  by  Technigaz  to  compute  the 
shear  and  bending  moment  curves  had  a  length  equal  to  the 
ship's  length-between-perpendiculars  which  was  275  meters. 
The  wave  height  was  equal  to  0,03  of  the  wave  length  or  8.4 
meters.   The  depth  of  water  at  frame  228  was  then  calculated 
by  adding  the  wave  height  to  the  Stillwater  draft.   The 
Smith  effect  was  also  included  in  these  calculations. 

An  interesting  load  curve  results  when  a  ship  carries 
cargo  in  large  spherical  tanks  as  with  this  ship.   The 
spherical  tank  is  supported  at  various  points  around  its 
equator,  thus  the  loads  are  transmitted  to  the  hull  of  the 
ship  in  a  circle.   The  load  curve  would  then  have  extremely 
high  loads  per  unit  length  in  the  areas  where  the  tanks  are 
adjoining  and,  fairly  low  loads  in  the  area  of  the  trans- 
verse centerlines  of  the  tanks.   Thus,  for  this  particular 
ship  the  load  curve  would  have  high  spikes  at  frames  61, 
129,  195,  261,  327,  and  389.   This  loading  curve  cannot  be 
accurately  approximated  by  a  series  of  straight  lines  indi- 
cating a  uniform  load  over  a  certain  segment  of  the  vessel's 
length  as  is  the  normal  design  practice.   The  load  curve  for 
this  particular  type  of  vessel  must  be  either  accurately 
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plotted  and  a  graphical  solution  obtained  or  some  type  of 
curve -fitting  computer  solution  performed. 

The  shear  and  moment  curves  provided  by  Technigaz  were 
calculated  using  the  dynamic  weight  of  the  cargo.   The 
dynamic  loads  in  this  region  were  calculated  by  determining 
the  vertical  accelerations  at  tank  3»  The  vertical  acceler- 
ation included  the  acceleration  due  to  gravity,  heave  accel- 
erations and  pitch  accelerations,  the  total  of  which  was 
13.21  meters/second.   These  accelerations  were  determined 
from  Series  60  data  that  gave  results  indicating  a  heave  of 
±7  meters  in  10.0  seconds  and  a  pitch  of  ±5»7  degrees  in 
10.0  seconds.   The  dynamic  force  was  then  calculated  by 
multiplying  the  mass  of  the  tank  and  LNG  times  the  ratio  of 
the  total  acceleration  to  the  acceleration  of  gravity.   This 
ratio,  which  was  the  number  of  "g"  forces  felt  at  tank  3»was 
1.4.   The  resulting  dynamic  load  due  to  the  tank  and  cargo 
was  20,200  metric  tons. 

Information  supplied  by  Technigaz  indicated  that  the 
steel  weight,  excluding  the  spherical  tank,  was  6,472  metric 
tons  between  frames  195  &nd  26l,  the  area  of  interest.   A 
quarter  tank  section  then  had  a  steel  weight  of  1,618  metric 
tons. 

The  only  additional  forces  that  had  to  be  applied  to 
the  model  were  the  boundary  forces  which  will  simulate  the 
remainder  of  the  ship.   As  is  indicated  by  Figure  19.  the 
after  portion  of  the  ship  may  be  simulated  by  the  applica- 
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tion  of  a  bending  moment  of  484,9^0  meter-tons  and  a  shear 
force  of  1,221  tons  at  frame  195»  the  free  end.   In  addition 
it  was  necessary  to  apply  a  shear  force  at  frame  228,  the 
transverse  plane  of  symmetry,  of  such  sign  and  magnitude  as 
to  maintain  the  model  in  rigid  body  equilibrium.   Or,  in 
other  words,  it  is  necessary  that  the  sum  of  the  forces  in 
the  Z  direction  be  equal  to  zero.   An  alternative  method 
used  by  some  investigators  in  this  field  to  insure  rigid 
body  equilibrium  is  to  fix  the  support  joints  in  the  trans- 
verse plane  of  symmetry  in  the  Z  direction.   This  method  is 
not  satisfactory  in  this  particular  case  since  knowledge  of 
deflections  of  the  members  is  the  primary  reason  for  the 
analysis  and  for  this  particular  ship  there  are  large  Z 
direction  deflections  in  the  plane  of  symmetry. 

First  Run — Loads 

The  beam  model  was  to  be  subjected  to  basically  four 
different  types  of  forces.   These  were  boundary  conditions, 
loads  due  to  water  acting  on  the  hull,  loads  due  to  the 
cargo,  and  loads  due  to  the  steel  hull  weight.   The  question 
was,  how  should  the  various  forces  be  distributed  to  accu- 
rately represent  the  actual  case. 

The  actual  steel  weight  of  the  quarter  section  of  the 
model  was  distributed  equally  among  the  joints  with  a  joint 
load  of  7.78  metric  tons  in  the  negative  Z  direction.   This 
was  considered  a  valid  distribution  since  there  were  a 
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greater  number  of  joints  in  locations  where  there  was  a 
large  amount  of  structural  steel.   The  concept  of  a  uniform 
force  in  the  negative  Z  direction  applied  to  each  member  was 
completely  infeasible  from  the  standpoint  of  additional  pro- 
gramming and  computer  time  required,  although,  this  would 
have  been  the  ideal  solution. 

The  dynamic  load  of  the  cargo  was  20,200  metric  tons  of 
5,050  tons  per  quarter  tank  section.   This  load  was  applied 
as  a  joint  load  at  each  of  the  vertical  members  of  the  sup- 
port platform  (third  deck  of  the  model).   The  load  was 
equally  distributed  among  the  joints  resulting  in  a  load  of 
505.0  tons  at  all  the  joints  except  joints  22  and  338  which 
were  in  the  planes  of  symmetry  and  consequently  had  one  half 
the  applied  load  of  the  other  joints  or  252.5  tons.   Sample 
Calculation  3  illustrates  the  method  of  calculation  of  the 
applied  loads. 

The  forces  due  to  the  water  acting  on  the  bottom  and 
side  of  the  hull  were  calculated  as  if  the  depth  of  the 
water  remained  constant  over  the  entire  quarter  tank  sec- 
tion,  This  assumption  was  slightly  in  error  in  that  the 
depth  was  slightly  less  at  the  ends  of  the  model  due  to  the 
fact  that  the  wave  crest  was  near  the  center  of  the  model. 
However,  it  was  felt  that  this  effect  was  negligible  due  to 
the  trocoidal  nature  of  the  wave  crest.   The  forces  on  the 
side  of  the  ship  were  idealized  as  linearly  varying  loads  on 
the  vertical  members.   The  bottom  forces  were  applied  to  the 
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model  as  uniform  forces — half  of  the  total  bottom  forces 
were  applied  to  the  longitudinal  members  and  half  to  the 
transverse  members.   A  ratio  of  ship  dimension  to  model 
dimension  was  used  in  the  calculation  of  these  forces  so 
that  the  force  actually  experienced  by  the  ship  is  applied 
to  the  model. 

The  remaining  loads  to  be  applied  were  the  boundary 
conditions.   The  bending  moment  of  484,9^0  meter-mtons  at 
frame  195  was  simulated  by  coupled  forces.   The  neutral  axis 
of  the  model  at  this  point  was  calculated  to  be  about  10.5 
meters  from  the  double  bottom.   Above  the  neutral  axis  a 
uniform  load  was  applied  in  the  positive  X  direction  to  the 
main  deck  between  the  inner  and  outer  web  frames.   In  addi- 
tion a  linear  load  was  applied  to  the  inner  and  outer  web 
bulkheads  from  the  main  deck  to  the  neutral  axis,  with  the 
third  and  fourth  decks  carrying  a  uniform  load  equal  to  the 
value  of  linear  distribution  at  the  intersection  of  the  webs 
and  deck  members.   The  other  half  of  the  couple  was  applied 
below  the  neutral  axis  in  the  negative  X  direction.   This 
was  accomplished  by  imposing  a  uniform  load  across  the 
entire  half  width  of  the  double  bottom,  with  a  linearly  de- 
creasing load  on  the  webs,  starting  at  the  double  bottom  up 
to  the  neutral  axis.   The  second  deck  carried  a  uniform  load 
equal  to  the  linear  load  applied  to  the  web  at  its  joints. 
See  Figure  20  for  a  graphical  representation  of  the  distri- 
bution. 
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F.  =  4.042  mtons/cnj 
F,  =  4.063  mtons/cm 
q     =   6l .05  nitons 


FIGURE   20 

FRAME  195  BOUNDARY  CONDITIONS 

HOGGING  CONDITION—RUN  I 
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No  forces  were  applied  to  members  nearer  the  center- 
line  than  the  inner  web  with  the  exception  of  the  double 
bottom.   The  reason  for  this  was  that  structural  members 
do  not  contribute  to  the  longitudinal  strength  of  the  vessel 
due  to  their  discontinuous  nature  in  way  of  the  holds  for 

i 

the  spherical  tanks.   The  inner  web  was  the  structural  sec- 
tion that  was  closest  to  the  centerline  of  the  vessel  and 
maintained  its  longitudinal  continuity  although  it  did  form 
part  of  the  tank  hold. 

The  shear  load  at  frame  195  was  1,221  mtons  or  610.5 
mtons  for  the  quarter  tank  section.   This  was  equally  dis- 
tributed as  a  joint  load  of  61.05  mtons  to  each  of  the  10 
joints  that  lie  in  the  intersections  of  the  plane  of  frame 
195  and  the  inner  and  outer  bulkheads  of  the  wing  tanks.   No 
shear  force  was  applied  to  the  vertical  bulkheads  near  the 
centerline  because  of  their  large  lightening  holes  and  poor 
vertical  strength  when  compared  to  the  inner  and  outer  web 
bulkheads. 

The  remaining  boundary  condition  is  the  shear  force  at 
frame  228.   This  force  was  determined  by  summing  the  verti- 
cal forces  applied  to  the  model.   The  negative  of  this 
summation  then  gave  the  force  necessary  to  keep  the  model 
in  equilibrium.   For  this  particular  loading  condition  this 
force  was  a  negative  1,389  mtons  which  was  distributed  as 
joint  loads  at  nodes  in  the  inner  and  outer  web  bulkheads 
at  frame  288.   See  Figure  21. 
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FIGURE  21 
FRAME  228  BOUNDARY  CONDITIONS 
HOGGING  CONDITION— RUN  I 
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First  Run — Results 

The  first  attempts  to  get  results  from  this  model  were 
unsuccessful  due  to  the  large  amount  of  core  and  facilities 
required  by  the  model.   Some  of  the  subprograms  that  make  up 
STRUDL  have  job  control  statements  inherent  in  them  that 
allow  for  a  certain  number  of  disks  to  be  in  operation  when 
that  subprogram  is  being  used.   This  particular  model  re- 
quired more  discs  than  were  provided  by  the  subprogram. 
Thus  it  was  necessary  to  enter  input  statements  which  over- 
rode some  of  the  job  control  statements  of  some  of  the  sub- 
programs and  allowed  the  computer  to  utilize  sufficient 
facilities. 

Results  for  the  first  run  were  finally  obtained  and  the 
deflections  appeared  satisfactory  when  looked  at  on  a  macro- 
scale.   However,  a  closer  examination  of  the  section 
stresses  indicated  several  problem  areas.   A  major  problem 
occurred  at  frame  195^   The  deflections  of  the  fourth  deck 
at  the  inner  web  in  the  positive  X  direction  were  large  over 
a  beam  length  as  indicated  by  Figure  22.   This  resulted  in 
all  members  eminating  from  joint  408  to  have  extremely  high 
bending  stresses.   Furthermore,  this  large  deflection  in  the 
X  direction  was  transmitted  to  the  circular  hold  bulkhead  at 
joint  283  causing  the  circular  support  platform  to  elongate 
with  resulting  high  stresses  in  all  members  in  the  vicinity. 

The  elongation  of  the  circular  hold  acting  in  conjunc- 
tion with  the  water  forces  on  the  side  of  the  ship,  caused  a 
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First  number  =  deflection  in  X  direction 
Second  number  =  deflection  in  Y  direction 
Third  number  =  deflection  in  Z  direction 


FIGURE  22 
DEFLECTIONS  OF  FRAME  195 
HOGGING  CONDITION— RUN  I 
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deflection  in  the  negative  Y  direction  of  the  circular 
support  platform  bulkhead  at  the  narrow  area  of  the  ship's 
hull  near  the  transverse  centerline  of  the  tank  (frames  216 
to  228) 0   This  effect  was  most  pronounced  at  the  second  deck 
level  and  caused  high  bending  stresses  in  the  vertical  mem- 
bers in  this  area. 

A  second  major  problem  area  was  the  main  deck,  again 
at  the  narrower  sections  between  frames  219  and  228.   High 
stresses  occurred  in  the  horizontal  members  due  to  the  high 
axial  stresses  combined  with  the  high  bending  stresses. 

The  third  major  problem  area  occurred  in  the  double 
bottom  at  the  junction  of  the  double  bottom  and  the  vertical 
strength  members  of  the  support  platform  again  in  the  vicin- 
ity of  frame  228.   The  joints  at  this  junction  had  large 
negative  Z  deflections,  which  caused  large  bending  stresses. 

A  fourth  problem  area  was  in  the  vicinity  of  the  origin 
in  the  double  bottom.   Members  eminating  from  the  origin  had 
a  positive  Z  deflection  for  a  radius  of  about  6  meters  in- 
stead of  the  expected  negative  deflection. 

Figures  23  and  24  show  the  longitudinal  stress  distri- 
bution at  frames  195  and  228  respectively.   Approximately 
2  5  per  cent  of  the  members  of  the  model  had  total  normal 
stresses  that  exceeded  the  yield  stress  of  mild  steel.   At 
this  point  a  complete  review  of  the  model  and  all  calcula- 
tions was  undertaken. 

The  first  discovery  was  that  an  error  of  about  5  per 
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FIGURE  23 
LONGITUDINAL  STRESSES—FRAME  195 
HOGGING  CONDITION— RUN  I 
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FIGURE  2k 
LONGITUDINAL  STRESSES— FRAME  228 
HOGGING  CONDITION— RUN  I 
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cent  had  been  made  in  the  calculation  of  the  bottom  water 
pressure  forces  acting  on  the  longitudinal  members.   This 
resulted  in  the  fact  that  the  sum  of  the  forces  in  the  Z 
direction  was  not  zero.   This  imbalance  in  forces  caused  a 
small  amount  of  rigid  body  motion  which  would  then  have 
caused  the  unexpected  deflections  in  the  vicinity  of  the 
origin  outlined  above  as  the  fourth  problem  area. 

Conversation  with  representatives  of  American  Technigaz 
brought  to  light  the  fact  that  the  dynamic  weight  of  steel 
had  been  used  in  the  calculation  of  the  shear  and  bending 
moment  curves  instead  of  the  static  weight.  This  factor  had 
not  been  taken  into  account  for  the  first  run.   This  cor- 
rection would  increase  the  negative  Z  direction  joint  loads 
that  simulated  the  steel  weight.  A  summation  of  the  forces 
in  the  Z  direction  then  indicated  that  the  shear  force  re- 
quired at  frame  228  to  prevent  rigid  body  motion  could  then 
be  reduced  from  1,389  mtons  to  7^2  mtons.  Consequently, 
correction  of  this  error  would  help  solve  the  problem  of 
high  stresses  in  the  double  bottom  at  the  intersection  of 
the  double  bottom  and  the  vertical  members  of  the  circular 
support  platform  which  was  described  earlier  as  problem 
three. 

It  was  obvious  that  the  extreme  elongation  of  the  cir- 
cular hold  and  the  inconsistent  X  direction  deflections  in 
frame  195  were  caused  by  the  boundary  forces  that  were  used 
to  simulate  the  bending  moment.   Thus  it  was  necessary  to 
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determine  whether  the  distribution  of  forces  that  was  used 
was  correct  and  if  not,  what  was  the  correct  distribution. 
Obvious  methods,  such  as  building  a  plastic  scale  model  of 
two  or  three  sections  of  the  ship,  loading  the  model  and 
getting  the  longitudinal  stress  distribution  at  frame  195i 
were  considered.   Another  method  that  could  have  been  used 
to  obtain  the  desired  information  was  a  large  macro-mesh 
model  consisting  of  plate  elements  which  could  be  developed 
using  the  finite  element  capability  of  STRUDL.   These  two 
methods  of  solution  had  to  be  discarded  because  of  time 
considerations.  A  literature  search  was  then  conducted  to 
try  to  find  longitudinal  stress  data  from  experiments  using 
plates  with  circular  cutouts. 

Papers  were  found  in  which  experiments  had  been  per- 
formed on  plates  subjected  to  a  uniform  tension  force  at  its 
ends  with  a  series  of  circular  cutouts  in  which  the  cutouts 
had  been  reinforced  with  a  combing.   This  would  have  approx- 
imated very  closely  the  main  deck  of  this  ship.   However, 
the  papers  that  contained  this  work  only  presented  stress 
results  for  a  transverse  section  through  the  center  of  the 
circle.   This  is  understandable  since  the  highest  stresses 
would  occur  at  this  point.   No  information  regarding  the 
stress  distribution  in  other  parts  of  the  plate  was  pro- 
vided. 

The  only  information  that  was  readily  available  was 
Kirch's  solution  for  a  central  circular  hole  in  an  infi- 
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nitely wide  plate  subjected  to  a  pure  tension  stress  (Ref . 
5).   The  results  of  his  solution  werei 

2  2       4 

o  =  %   (1  -  T-%   ) ■ .+  § ■  (1  -  4  T~%    +  3  t;  )  cos  29 
r   2      r2     2        r2      r 


2  4 

ae  =  I  (1  +  ~l  }  "  I  (1  +  3  ~4  }  cos  2e 


where    a  =  tension  stresses  to  which  plate  is  subjected 
a     =  radial  stresses 

aa   =  tangential  stresses 

r  =  radius  of  hole 
o 

r  =  radius  at  which  calculation  is  being  made 

9  =  angle  r  makes  with  longitudinal  line  through 

center  of  the  circle 

Subsequent  to  the  calculation  of  the  radial  and  tangential 
stresses  through  a  transverse  section  20. 34  meters  ahead  of 
the  center  of  the  circle  it  was  possible  to  calculate  the 
distribution  of  the  longitudinal  stresses  by  the  equation: 

a=o  cose  +  0asin6 
x    r         9 

where    o  =  longitudinal  stress 

A. 

The  result  of  the  calculations  utilizing  the  above 
equations  is  shown  in  Figure  25  with  the  stresses  indicated 
as  a  fraction  of  the  stress  at  what  would  be  the  edge  of  the 
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FIGURE  2  5 
LONGITUDINAL  STRESS  DISTRIBUTION  IN  MAIN  DECK  CALCULATED 
ACCORDING  TO  KIRSCH'S  SOLUTION  FOR  AN  INFINITE  PLATE 
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main  deck. 

There  are  several  important  differences  between  the 
actual  main  deck  plating  and  the  idealized  plate  with  the 
circular  cutout.   First,  in  the  actual  situation  the  main 
deck  is  not  infinite  in  the  horizontal  plane  as  was  the 
above  solution.   Secondly,  Kirsch's  solution  did  not  take 
into  account  the  fact  that  the  actual  case  had  an  inner  web 
bulkhead  running  longitudinally  or  that  it  had  a  vertical 
wall  around  the  circumference  of  the  hole.   In  addition  the 
solution  for  the  infinite  plate  with  the  single  circular 
cutout  did  not  take  account  of  any  change  in  the  stress 
distribution  that  a  series  of  holes  in  the  longitudinal 
direction  would  cause. 

In  spite  of  the  many  differences  between  the  actual 
main  deck  and  the  plate  with  the  circular  cutout  and  with 
no  way  of  quantitatively  determining  the  effects  of  these 
changes,  it  was  felt  that  the  calculated  distribution  would 
be  helpful  in  determining  the  correct  boundary  conditions 
that  should  be  imposed  on  the  beam  model. 

It  was  readily  apparent  from  a  comparison  of  Figures 
23  and  25  that,  in  no  way,  did  the  distribution  of  longitu- 
dinal stresses  that  was  the  result  from  the  first  run  com- 
pare with  the  distribution  obtained  from  Kirsch's  solution. 
For  example,  the  longitudinal  stress  in  the  plate  calcula- 
tion at  the  inner  web  bulkhead  was  19  per  cent  of  the  load 
at  the  point  corresponding  to  the  side  of  the  ship  while  the 
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results  from  the  first  run  indicated  that  the  stress  at  the 
inner  web  was  88  per  cent  of  that  at  the  main  deck  edge. 

Thus  it  was  obvious  that  a  more  realistic  application 
of  the  bending  moment  boundary  conditions  was  necessary  to 
obtain  any  meaningful  results.   This  new  distribution  of 
forces  would  have  to  put  a  greater  load  on  the  part  of  the 
main  deck  that  is  longitudinally  continuous  as  well  as  a 
greater  share  on  the  outer  hull.   Such  a  refined  approach 
would  tend  to  reduce  the  problems  encountered  on  the  fourth 
deck  at  the  intersection  of  the  circular  hold  and  the  inner 
web  which  resulted  in  an  extreme  elongation  of  the  circular 
hold.  A  more  realistic  distribution  would  also  tend  to  re- 
duce the  high  stresses  that  occurred  at  the  narrow  portion 
of  the  main  deck  at  frame  228. 

Second  Run — Loads 

The  decision  was  made  to  run  the  STRUDL  program  a  sec- 
ond time  for  the  hogging  condition  making  the  changes  to  the 
loads  that  were  necessary.   The  first  correction  made  was  to 
change  the  joint  loads  due  to  the  ships  steel  weight  from 
7.78  metric  tons/joint  to  10. 89  metric  tons/joint.   This  new 
weight  was  the  result  of  using  a  "g"  force  of  1.4  as  was 
previously  discussed.   This  change  also  made  necessary  a 
change  in  the  shear  force  at  frame  228.   A  summation  of 
forces  in  the  Z  direction  indicated  that  the  shear  force  per 
joint  should  be  reduced  to  7^.2  mtons/ joint  as  indicated  in 
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Figure  26.   The  force  due  to  the  water  acting  on  the  long- 
itudinal members  of  the  double  bottom  was  also  corrected  for 
this  run. 

The  largest  problem  for  the  second  run  was  to  determine 
the  distribution  of  longitudinal  forces  that  would  accurate- 
ly simulate  the  bending  moment  and  give  a  reasonable  distri- 
bution of  longitudinal  stresses  in  the  main  deck.   After  a 
careful  reanalysis  it  was  decided  that  no  forces  should  be 
applied  to  the  horizontal  members  making  up  the  second, 
third,  or  fourth  deck.   The  reason  for  this  was  that  these 
decks  had  almost  no  longitudinal  strength  due  to  the  large 
lightening  holes  in  them.   The  elimination  of  forces  on 
these  members  would  be  a  conservative  estimate  when  the  ship 
is  looked  at  as  a  whole,  since  the  decks  would  assume  a 
small  portion  of  the  longitudinal  load.   In  order  to  try  to 
obtain  a  longitudinal  stress  distribution  in  the  main  deck 
a  linearly  varying  load  was  applied  to  the  transverse  mem- 
bers of  the  main  deck  of  frame  195  between  the  inner  and 
outer  web  frames.   The  value  of  the  force  at  the  inner  web 
was  half  of  that  at  the  edge  of  the  ship  as  indicated  in 
Figure  27.   In  addition,  a  linear  load  was  applied  to  the 
vertical  member  of  the  inner  and  outer  web  bulkheads,  with 
the  maximum  value  being  the  force  on  the  main  deck  or  the 
double  bottom  and  a  zero  force  at  the  neutral  axis. 

Second  Run — Results 

The  results  from  the  second  run  seemed  much  more  rea- 
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FIGURE  26 
FRAME  228  BOUNDARY  CONDITIONS 
HOGGING  CASE—RUN  II 
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FIGURE  27 
FRAME  195  BOUNDARY  CONDITIONS 
HOGGING  CASE— RUN  II 
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sonable  than  had  those  of  the  first  run.   The  positive  Z 
direction  deflection  of  the  double  bottom  near  the  origin 
had  completely  disappeared  with  the  elimination  of  the  rigid 
body  motion  that  had  been  caused  by  the  miscalculation  of 
the  longitudinal  forces.   However,  as  can  be  seen  from 
Figure  28  the  longitudinal  stress  distribution  on  the  main 
deck  at  frame  195  still  did  not  come  close  to  the  distribu- 
tion obtained  for  the  exact  solution  of  the  infinite  plate 
with  the  circular  hole.   The  results  from  Run  II  indicate 
that  the  stress  at  the  inner  web  is  approximately  54  per 
cent  of  the  longitudinal  stress  at  the  edge  of  the  main 
deck.   This  5^  per  cent  represented  an  improvement  over  the 
88  per  cent  obtained  in  the  first  run  but  still  falls  short 
of  the  desired  result  of  19  per  cent  obtained  from  the  in- 
finite plate  solution.   It  was  obvious  that  the  forces 
applied  to  the  main  deck  that  are  part  of  the  couple  used  to 
simulate  the  bending  had  to  be  redistributed  again. 

The  problem  on  the  fourth  deck  where  the  inner  web 
joined  the  circular  hold  wall  was  greatly  improved,  however, 
the  stress  in  some  members  was  still  excessive.   The  bending 
stresses  remained  the  reason  for  this  problem  with  members 
42?,  ^28,  and  429  being  pulled  out  of  the  circular  shape  by 
member  409.   The  poor  distribution  of  the  bending  moment 
couple  forces  that  was  discussed  in  the  previous  paragraph 
would  be  one  of  the  primary  reasons. 

As  can  be  seen  from  Figure  29  the  longitudinal  stresses 
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FIGURE  28 
LONGITUDINAL  STRESSES  AT  FRAME  195 
HOGGING  CONDITION— RUN  II 
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FIGURE  29 
LONGITUDINAL  STRESSES— FRAME  228 
HOGGING  CONDITION— RUM  II 
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in  this  run  still  exceeded  the  yield  stress  in  the  main  deck 
at  frame  228.   This  problem  was  discussed  at  a  conference 
with  Technigaz  representatives.   At  this  conference  it  was 
discovered  that  there  had  been  a  communication  problem  con- 
cerning a  revision  to  the  blueprints  from  which  this  work 
was  done.   The  revision  concerned  an  additional  longitudinal 
bulkhead  under  the  main  deck  that  had  been  added  for  struc- 
tural purposes.   Runs  I  and  II  were  made  using  20  milli- 
meters as  the  thickness  of  the  bulkhead.   The  correct  thick- 
ness of  this  bulkhead  was  ^5  millimeters.   Thus  it  was  nec- 
essary to  make  a  third  run  to  obtain  valid  results. 

Third  Run — Loads 

Most  of  the  loads  for  this  run  remained  the  same  as  for 
the  second  run,  the  only  exception  being  the  forces  making 
up  the  upper  half  of  the  couple  that  was  simulating  the 
bending  moment  at  frame  195*   The  distribution  of  these 
forces  was  changed  such  that  the  linearly  varying  force  be- 
tween the  inner  and  outer  web  had  a  value  at  the  inner  web 
that  was  one  third  of  the  value  at  the  deck  edge.   The  ver- 
tical members  in  the  inner  and  outer  webs  were  also  subjec- 
ted to  an  X  direction  linearly  varying  load  that  had  its 
maximum  value  at  the  main  deck  and  had  a  value  of  zero  at 
the  neutral  axis.   The  boundary  conditions  at  frame  195  are 
shown  graphically  in  Figure  30.   The  boundary  conditions  at 
frame  228  remain  the  same  as  shown  in  Figure  26. 


-72- 


Fb  = 

q  = 


7.799  mtons/cm 
4,156  mtons/cm 
61.05  mtons 


FIGURE  30 
FRAME  195  BOUNDARY  CONDITIONS 
HOGGING  CASE— RUN  III 
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In  addition,  the  properties  of  members  501  through  518 
were  recalculated  to  account  for  the  error  in  the  bulkhead 

thickness. 

Third  Run — Results 

The  results  from  the  third  run  were  better  than  the 
first  or  second  runs.   The  shape  of  the  longitudinal  stress 
distribution  in  the  main  deck  at  frame  195  (Fig.  31)  con- 
formed more  closely  to  that  of  the  infinite  plate  solution. 
The  value  of  the  longitudinal  stress  at  the  inner  web  was 
4  5  per  cent  of  that  at  the  deck  edge.  Although  this  is  not 
the  value  obtained  from  the  infinite  plate  model,  time  con- 
siderations as  well  as  the  high  cost  of  each  run  indicated 
that  further  experimentation  with  the  boundary  conditions 
simulating  the  bending  moment  on  the  upper  half  of  the  ship 
was  unwarranted.   This  decision  was  based  on  two  considera- 
tions.  The  first  was  the  uncertainty  of  the  actual  longi- 
tudinal stress  distribution  -  the  infinite  plate  had  many 
drawbacks  as  was  previously  cited.   The  second  consideration 
was  that  the  distribution  of  the  forces  on  the  main  deck 
would  have  very  little  effect  on  the  vertical  deflections  of 
the  support  platform  -  the  primary  desired  result.   Figures 
33  through  37  show  the  deflections  obtained  for  the  third 
run. 

The  maximum  vertical  deflection  of  any  point  in  re- 
lation to  any  other  point  on  the  circular  platform  for  the 
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FIGURE  31 
LONGITUDINAL  STRESSES—FRAME  195 
HOGGING  CONDITION—RUN  III 
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FIGURE  32 
LONGITUDINAL  STRESSES  —FRAME  228 

HOGGING  CONDITION—RUN  III 
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FIGURE   33 
BOTTOM  DECK   DEFLECTIONS 
HOGGING   CONDITION— RUN   III 
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First  number  =  deflection  in  X  direction 
Second  number  =  deflection  in  Y  direction 
Third  number  =  deflection  in  Z  direction 
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FIGURE  Jk 

SECOND  DECK  DEFLECTIONS 

HOGGING  CASE— RUN  III 
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THIRD  DECK  DEFLECTIONS 
HOGGING  CASE— RUN  III 
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First  number  =  deflection  in  X  direction 
Second  number  =  deflection  in  Y  direction 
Third  number  =  deflection  in  Z  direction 
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FIGURE  36 
FOURTH  DECK  DEFLECTIONS 

HOGGING  CONDITION— RUN  III 
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First  number  =  deflection  in  X  direction 
Second  number  =  deflection  in  Y  direction 
Third  number  =  deflection  in  Z  direction 
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FIGURE  37 
MAIN  DECK  DEFLECTIONS 
HOGGING  CONDITION— RUN  III 
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first,  second,  and  third  runs  were  1.0  cm,  0.9  cm,  and  0.6 
cm,  respectively.   Even  with  major  changes  in  boundary  condi- 
tions (between  run  I  and  run  II)  and  changes  in  the  longitu- 
dinal strength  of  the  main  deck  (between  run  II  and  run  III) 
this  deflection  difference  remained  fairly  constant. 

Further  support  for  the  argument  that  the  distribution 
of  forces  would  have  little  effect  on  the  vertical  deflec- 
tions can  be  seen  by  looking  at  the  magnitude  of  the  verti- 
cal deflections  between  runs  II  and  III  in  which  the  verti- 
cal loads  were  constant.   The  smallest  vertical  deflection 
on  the  support  platform  occurred  at  joint  22  and  was  -6.5  cm 
and  -6.6  for  runs  I  and  II  respectively,  while  the  largest 
occurred  at  joint  338  and  was  <-7.4  cm  and  -7.2  cm  respec- 
tively.  The  conclusions  being  that  changes  in  the  distri- 
bution of  forces  simulating  the  bending  moment  and  even 
changes  in  the  longitudinal  strength  of  the  main  deck  have 
little  effect  on  the  vertical  deflections  of  the  support 
platform. 

The  same  is,  naturally,  not  true  for  the  transverse  or 
longitudinal  deflections  since  the  boundary  conditions  will 
greatly  affect  the  longitudinal  elongation  and  subsequent 
transverse  contraction  of  the  circular  support  platform  when 
the  ship  is  in  a  hogging  situation.   The  maximum  deflection 
in  the  X  direction  of  the  support  joints  occurred  in  the 
longitudinal  direction  and  was  due  to  the  elongation  of  the 
circular  support  platform.   This  joint  deflection  was  4.4 
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cm,  3.3  cm,  and  2.5  cm  for  the  first,  second,  and  third 
runs  respectively. 

The  largest  Y  direction  deflection  occurred  at  the 
transverse  centerline  of  the  tank.   This  contraction  of  the 
circular  platform  in  the  transverse  direction  was  directly 
linked  to  the  elongation  in  the  longitudinal  direction.   The 
maximum  magnitude  of  this  transverse  deflection  was  4.1  cm, 
2.8  cm,  and  2.1  cm  for  the  three  runs. 

The  third  run  had  the  most  realsitic  longitudinal  dis- 
tribution.  Thus,  the  maximum  expected  deflection  of  any 
joint  normal  to  the  circular  shape  for  this  loading  condition 
was  2.5  cm.   It  was  felt  that  this  number  was  conservative 
and  could  be  used  for  further  work  on  the  preliminary  de- 
sign. 

The  areas  of  high  stresses  that  occurred  in  the  third 
run  were  in  basically  the  same  areas  as  in  the  two  preceed- 
ing  runs,  although  the  magnitude  of  the  stresses  were  lower. 
The  maximum  longitudinal  stress  occurred  at  the  narrow  por- 
tion of  the  main  deck  as  would  be  expected.   The  longitudi- 
nal stresses  when  combined  with  the  bending  stresses  re- 
sulted in  a  total  member  stress  that  was  greater  than  yield 
for  the  beams  in  this  area. 

Vertical  members  in  this  same  area  near  the  narrow  por- 
tion of  the  circular  hold  between  the  lower  two  decks  had 
combined  stresses  that  exceeded  the  yield  stress.   The 
primary  cause  of  this  was  the  bending  stresses  which  result- 
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ed  from  a  combination  of  loads  involving  the  transverse  con- 
traction  of  the  circular  hold  area  and  the  loads  due  to  the 
water  pressure  acting  on  the  side  of  the  hull.   In  addition 
vertical  members  of  the  circular  support  platform  at  the 
ship's  longitudinal  centerline  indicated  high  bending 
stresses  due  to  the  longitudinal  elongation  of  the  circular 
support  platform.   The  same  elongation  produces  high  bending 
stresses  in  members  427  and  428. 

It  was  felt  that  the  fact  that  beam  elements  were  being 
used  to  model  plates  in  the  outer  skin  of  the  ship  was  the 
cause  of  the  high  bending  stresses  in  this  region.   The 
linear  force  distribution  simulating  the  water  pressure 
caused  the  members  to  bow  in,  resulting  in  high  stresses. 
However,  in  the  actual  case  there  would  be  another  point  of 
support  in  the  middle  of  the  beam,  cutting  the  amount  of 
bowing  as  was  previously  explained  in  the  mathematical  model 
chapter. 

However,  it  was  the  opinion  of  the  writer  that  steps 
should  be  taken  to  insure  that  the  bending  stresses  due  to 
the  elongation  and  contraction  of  the  circular  hold  are  sub- 
stantially smaller  than  the  results  indicated  for  this 
model.   The  same  is  true  for  the  normal  stresses  in  the  main 
deck  at  the  narrow  portion  of  the  hull  near  the  transverse 
centerline  of  the  tank. 

The  increase  in  the  area  of  members  501  and  510  should 
have  brought  about  a  proportional  decrease  in  the  longitudi- 
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rial stress  level  of  the  main  deck  at  frame  228.  This  stress 
reduction  could  not  be  checked  directly  for  each  member  be- 
cause the  distribution  of  boundary  forces  in  the  main  deck 
at  frame  195  had  been  changed  between  run  II  and  run  III. 
However,  the  axial  force  on  each  member  was  one  of  the  out- 
puts of  this  program.   These  forces  were  used  as  correction 
factors  to  obtain  the  stress  level  of  the  members  if  the 
same  distribution  had  been  used  for  run  II  as  was  used  for 
run  III. 

This  was  accomplished  by  using t 

F  =  aA 

* 
and      F2  =  a2A2  F_  =  OgAg 


since  Ap  =  A2 


F2/a2  =  F3/a2 


°2  ■  a2  (W 


where    a2  =  actual  longitudinal  stress  obtained  for  run  II 

0?  =  longitudinal  stress  that  would  have  been  ob- 
tained for  run  II  if  the  boundary  conditions 
of  run  III  had  been  used 

F2  =  axial  force  in  member  for  run  III 
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F-  =  axial  force  in  member  for  run  III 

A„  =  cross  sectional  area  of  member  for  run  II 


Since  the  stresses  and  area  can  be  considered  inversely 
proportional  if  a  constant  force  is  applied,  the  graph  of 
the  stress  versus  area  of  a  member  would  be  approximately 
linear.   Thus,  the  information  from  runs  II  and  III  can  be 
used  to  determine  the  area  required  for  any  desired  stress 
level  when  the  vessel  is  in  the  hogging  condition.   See 
Figure  38. 
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FIGURE   38 

RELATIONSHIP  OF  AREA  AND  STRESS  LEVEL 

IN  MEMBERS  501  AND  510  FOR  HOGGING  CASE 
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THE  ANALYSIS— SAGGING  CONDITION 

General  Discussion 

The  last  case  to  be  investigated  was  with  the  ship  in 
an  upright  position,  fully  loaded  with  the  wave  crests  at 
the  end  of  the  ship  and  the  wave  hollow  at  amidships.  Again 
the  wave  height  was  taken  to  be  0.03  of  the  wave  length  or 
8.4  meters.   The  draft  at  frame  228  was  then  calculated  by 
subtracting  the  wave  height  from  the  stillwater  draft,  this 
resulted  in  a  draft  of  6.80  meters. 

The  shear  and  moment  curve  for  this  loading  condition 
were  also  porvided  by  Technigaz.  The  bending  moment  at 
frame  195  was  -309,800  meter-tons  while  the  shear  load  was 
296  mtons  as  indicated  by  Figure  39.  Since  the  dynamic 
loads  would  be  slightly  less  than  the  static  loads  for  the 
sagging  condition,  the  static  loads  were  used  to  give  con- 
servative results. 

First  Run — Loads 

Again  the  beam  model  was  subjected  to  boundary  loads, 
loads  due  to  water  acting  on  the  hull,  cargo  loads  and  loads 
due  to  the  steel  weight.   The  static  steel  weight  of  7.78 
mtons/joint  was  applied  to  each  joint  of  the  model. 

The  cargo  load  was  14,428  metric  tons  or  3t607  metric 
tons  per  quarter  tank  section.   The  load  was  equally  distri- 
buted among  the  joints  resulting  in  a  load  of  360.7  metric 
tons  at  all  the  joints  except  joints  22  and  338  which  were 
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in  the  planes  of  symmetry  and  had  one  half  the  applied  load 
of  the  outer  joints  or  180.4  metric  tons. 

The  forces  due  to  the  water  acting  on  the  bottom  and 
side  hull  were  calculated  as  if  the  depth  of  water  remained 
constant  over  the  entire  quarter  tank  section.   The  forces 
on  the  side  of  the  ship  were  idealized  as  linearly  varying 
loads  on  the  vertical  members.   As  in  the  hogging  case  the 
bottom  forces  were  applied  to  the  model  as  uniform  forces — 
half  of  the  total  forces  to  the  longitudinal  members  and 
half  to  the  transverse.   See  Sample  Calculation  3. 

The  bendjng  moment  was  simulated  by  coupled  forces 
taking  the  neutral  axis  to  be  10. 5  meters  from  the  double 
bottom.   The  bending  moment  of  309,800  meter-mtons  was  dis- 
tributed as  it  was  for  Run  III  of  the  hogging  case.  Above 
the  neutral  axis  a  linear  load  was  applied  in  the  negative 
to  the  transverse  member  of  the  main  deck  between  the  inner 
and  outer  frames.   The  magnitude  of  the  force  at  the  inner 
web  was  one  third  of  the  force  at  the  deck  edge.   A  linear 
load  was  also  applied  to  the  inner  and  outer  bulkheads  from 
the  main  deck  to  the  neutral  axis.   The  other  half  of  the 
couple  was  applied  below  the  neutral  axis  in  the  positive  X 
direction.   A  uniform  load  was  imposed  across  the  entire 
half  width  of  the  double  bottom,  with  a  linearly  decreasing 
load  on  the  webs,  starting  at  the  double  bottom  up  to  the 
neutral  axis.   See  Figure  40  for  a  graphical  representation. 

The  shear  load  of  296  mtons  for  the  whole  cross  section 


-90- 


"t"              ,   mtons/cm 
7     =  2.655  mtons/cm 

q  =  14.7  mtons 

FIGURE  kO 

FRAME  195  BOUNDARY  LOADS 

SAGGING  CASE— RUN  I 
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or  148  mtons  for  the  quarter  tank  model  was  equally  distri- 
buted as  a  joint  load  of  14.8  mtons  acting  in  the  positive  Z 
direction  on  each  of  the  ten  joints  that  lie  in  the  inter- 
section of  frame  195  and  the  inner  and  outer  web  bulkheads. 
The  one  remaining  boundary  condition,  the  shear  at  frame 
228,  was  determined  by  summing  the  vertical  forces  applied 
to  the  model  and  then  changing  the  sign  of  this  summation 
thus  causing  the  model  to  be  in  equilibrium  in  the  vertical 
direction.   This  shear  at  frame  228  was  determined  to  be 
1598  mtons.   This  resulted  in  a  joint  load  of  159.8  mtons 
for  each  of  the  joints  in  the  inner  and  outer  web  bulkheads 
of  frame  228  (Fig.  41). 

First  Run — Results 

Generally,  the  deflections  and  stresses  for  the  sagging 
case  were  less  than  for  the  hogging  case.  As  would  be  ex- 
pected the  directions  of  the  longitudinal  and  transverse 
deflections  were  the  opposite  of  those  in  the  hogging  case. 
The  circular  tank  hold  contracted  in  the  longitudinal  direc- 
tion and  elongated  in  the  transverse  direction. 

The  shape  of  the  longitudinal  stress  distribution  in 
the  main  deck  conformed  almost  exactly  to  the  distribution 
obtained  in  the  third  run  of  the  hogging  case  although  the 
sign  was  different  because  the  deck  was  in  compression  for 
this  case  (Fig.  42).  The  longitudinal  stress  at  the  inner 
web  was  44  per  cent  of  the  stress  at  the  deck  edge.   This 
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FIGURE  kl 
FRAME  228  BOUNDARY  LOADS 
SAGGING  CASE—RUN  I 
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LONGITUDINAL  STRESSES—FRAME  195 
SAGGING  CASE— RUN  I 
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LONGITUDINAL  STRESSES— FRAME  228 
SAGGING  CASE— RUN  I 
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result  is  very  close  to  the  45  per  cent  obtained  for  the 
hogging  case.   If  the  main  deck  had  been  an  infinite  plate 
subjected  to  pure  compression  or  tension  the  distributions 
would  have  been  exactly  the  same. 

The  deflections  in  the  sagging  case  were  of  a  smaller 
magnitude  than  the  other  runs  as  can  be  seen  from  Figures 
44  through  48.   Figure  46  indicates  that  the  maximum  deflec- 
tion of  any  point  in  relation  to  any  other  point  on  the  cir- 
cular support  platform  was  0.3  cm.   The  maximum  deflection 
of  the  joints  in  the  support  platform  carrying  the  cargo 
load  was  0.7  cm  and  occurred  at  the  transverse  tank  center- 
line  due  to  the  transverse  elongation.   The  maximum  contrac- 
tion of  the  support  platform  occurred  near  the  ship's  longi- 
tudinal centerline  and  had  a  magnitude  of  0.5  cm. 

There  were  only  four  members  that  had  a  total  normal 
stress  ir  excess  of  the  yield  stress.   Two  of  these  members, 
609  and  627,  are  vertical  members  in  frame  lf5  and  the  high 
local  bending  stresses  are  due  to  the  bowing  caused  by  the 
application  of  the  boundary  conditions.   The  high  stresses 
are  completely  unrealistic  and  are  of  no  concern  in  the 
actual  ship.   The  other  two  members,  672  and  687,  have  high 
stresses  due  to  bending.   This  bending  was  caused  by  longi- 
tudinal contraction  of  the  circular  tank  hold  in  the  case  of 
member  672  and  the  transverse  elongation  in  the  case  of  mem- 
ber 687.   The  stresses  in  these  two  members  bears  further 
investigation.   However,  it  was  the  opinion  of  the  writer 
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First  number  =  deflections  in  X  dierection 
Second  number  =  deflections  in  Y  direction 
Third  mmber  =  deflections  in  Z  direction 
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FOURTH  DECK  DEFLECTIONS 

SAGGING  CASE— RUN  I 
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MAIN  DECK  DEFLECTIONS 

SAGGING  CASE— RUN  I 
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that  these  high  bending  stresses  are  the  result  of  the  beam 
model  and  are  not  realistic. 
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CONCLUSIONS 

i 

1.  The  ICES  STRUDL  program  can  be  effectively  utilized 
in  the  preliminary  design  phase  to  determine  stresses  and 
deflections  of  ships  of  unusual  form  for  which  traditional 
design  criteria  is  unsatisfactory. 

2.  The  beam  element  model  used  gave  good  results  for 
the  deflections  of  the  circular  support  platform  in  the 
vertical  direction.   These  results  were  basically  indepen- 
dent of  end  effects  due  to  the  application  of  loads  simu- 
lating the  remainder  of  the  ship.   The  maximum  vertical  de- 
flection of  any  point  in  relation  to  any  other  point  on  the 
circular  support  platform  was  0.6  centimeters  for  the  con- 
ditions investigated. 

3.  The  results  for  the  magnitude  of  the  transverse  and 
longitudinal  deflections  of  the  circular  platform  were  in- 
fluenced by  end  conditions,  particularly  the  distribution  of 
forces  simulating  the  bending  moment.   Upper  limits  were 
obtained  for  these  magnitudes  using  a  conservative  distri- 
bution of  forces. 

k.      The  longitudinal  stresses  were  also  influenced  by 
the  end  condition.  However,  it  was  possible  to  obtain  an 
upper  limit  on  the  maximum  value  using  a  conservative  dis- 
tribution of  boundary  conditions.   This  maximum  longitudinal 
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stress  occurred  in  the  main  deck  at  the  narrow  portion  of 
the  hull  as  would  be  expected.   The  stress  level  in  this 
area  was  of  sufficient  magnitude  to  warrant  consideration  of 
the  use  of  high  strength  steel  in  this  area. 

5.   The  local  member  bending  stress  results  obtained 
did  not  accurately  reflect  the  actual  case.   This  is  due  to 
the  inability  of  beam  elements  to  accurately  model  plates 
on  a  local  scale.   However,  indications  of  areas  with  high 
bending  stresses  were  obtained  and  these  should  be  subjected 
to  further  investigation.   Possible  reduction  of  cutout  sizes 
or  additional  flange  reinforcement  at  the  holes  are  possible 
cures  if  needed. 
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RECOMMENDATIONS 

i 

This  thesis  investigated  the  stresses  and  deflections 
for  the  hogging  and  sagging  condition  with  the  ship  in  an 
upright  mode.   It  is  recommended  that  additional  work  be 
carried  out  to  determine  the  ships  reaction  to  rolling  and 
to  torsion.   The  torsional  stresses  incurred  when  this  ship 
is  heading  obliquely  to  a  wave  would  be  appreciable  due  to 
the  small  amount  of  main  deck  area. 

A  further  recommendation  is  that  the  end  effects  be 
eliminated  in  any  future  work  with  this  model  or  any  other 
single  tank  section  model.   This  can  be  accomplished  by  any 
of  three  wayst 

(1)  The  model  could  be  expanded  in  the  longitudinal 
direction,  thus  eliminating  the  end  effects  on  the  quarter 
tank  section  of  interest. 

(2)  A  scale,  strain  guage  equipped  model  could  be  con- 
structed so  that  the  correct  distribution  of  longitudinal 
stresses  in  the  main  deck  and  web  bulkheads  at  frame  195  can 
be  determined. 

(3)  A  macro-scale  finite  element  program  could  be 
developed  that  would  give  as  output  the  longitudinal  stress 
distribution  at  the  main  deck. 
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APPENDIX  A 
Sample  Calculations 
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SAMPLE  CALCULATION  #1 


EFFECTIVE  THICKNESS  OF  PLATE 


Double  Bottom  -  Transverse  Cross  Section 


L    [     I 


r   r  r  \\  r  r  r 


3A8 


L    I    I 


-*- 


3.18 


T,=  /6  n»»n 


T.«  '8.*  »« 


stiffner  spacing  =  0.82  meters 

frame  spacing  =  3»28  meters 

cross  sectional  area  of  400  X  Ik   HP  stiffners  =  6125  mm' 

cross  sectional  area  of  400  X  15  HP  stiffners  =  8450  mm' 


At  -  area  of  stiffner  x  (n0#  stiffner  between  frames) 
frame  spacing 


effective 


'effective 


=  16  +  £|fo  X  3  =  21.6  mm 
=  18'5  +  §2§§  x  3  =  26,2  mm 


z 
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SAMPLE  CALCULATION  #2 


MEMBER  PROPERTIES 


Double  Bottom  -  Transverse  Cross  Section 


L,  s  1.29 


,.  T,?  18.*  mm 


-»  y 


lx'  *A8  m 


local  coordinates 


T,  * au  **v 

4 


t  t 


»•>= 


A  * 


3.S*  -m 


The  double  bottom  longitudinals  may  be  idealized  as 
beams  as  indicated  by  the  solid  lines  above • 

AX  =  cross-sectional  area 

=  (3.28)  (.0216)  +  (2.50)  (.0185)  +  (3.28)  (.0262) 
=  .2030  m2 

AY  =  shear  area  in  Y  direction  =  area  of  flanges  of  beam 
=  (3.28)  (.0216)  ♦  (3.28)  (.0262) 
=  .1567  m2 

AZ  =  shear  area  in  Z  direction  =  area  of  web  of  beam 
=  (2.50.)  (.0185) 
=  .0462  m2 
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IX  =  torsional  rigidity 
=  (1/3)   L(teff)3 

=  (1/3)   (3.28)(.02l6)3*  (3.28)(.0262)3+  (2.50) ( .0185)3 

=  .000036  m 

The  effective  breadth  of  the  plating  is  equal  to  60  times 

the  thickness  of  the  plate 

in  the  y  direction 

Li    =  60(.02l6)  =  1.296  m 
eff 

L2    =  60(.0262)  =  1.572  m 
eff 

in  the  z  direction 

Lo    =  60(.0185)  =  1.11  m 
Jeff 

It  is  necessary  to  determine  the  y  and  z  neutral  axes  in 
order  to  calculate  the  moments  of  inertia  and  section  mod- 
uli for  the  y  and  z  axes. 

y  neutral  axis  -  it  is  necessary  to  use  L<    and 

eff 

Lo    to  calculate  the  y  neutral  axis 
eff 
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■t,  •   M«4* 


Y*«! 


7  =  I8.t»* 


'<« 


Xt.L   -m» 


k,'-  J'f* 


1   Ta 


=  14.1*w 


f    «« 


YNA  =  distance  from  bottom  of  member  to  neutral  axis 
Y  =  distance  from  the  y  neutral  axis  to  the  extreme  section 

YNA  =  2  (areas) (moment  arms) 
L  areas 

_    (1.296H.0216H2.50)   +    (2.  50)(. 018SH2.  50/2) 
(1.296)(.02l6)   +    (2.50)(.0185)   ♦    (1.572)(.0262) 

=   1.148  m 


IY  = 


SY  = 


moment  of  inertia  about  y  axis  (using  L*  and  L?    ) 

Leff       eff 

(1/12)(.0185)(2.50)3  +  (.0216)(2.50  -  1.107)2 

+  (.0262)(l.572)(1.107)2 

k 

.1322  m 

section  modulus  about  y  axis 

II 
Y 

.1322 
1.352 


=    .0978    m- 
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u,  *  *.***• 


7J  '  I8.S  -* 


*  1   e  ll.Un 


L.    =  1,11   *» 


1   f     c  3fc,2  >»i 


* Lj*  ?.3»  *> » 


«     «** 


Z  =  distance  from  the  z  neutral  axis  to  the  extreme  section 

The  z  neutral  axis  is  a  vertical  line  through  the  horizontal 
centers  of  the  three  sections  of  the  member. 


IZ  =  moment  of  inertia  about  Z  axis  (use  L~    ) 

^eff 

=  (1/I2)(.02l6)(3.28)  +  (l/l2)(.0262)(3.28)-;, 


SZ  = 


.1406  m 
IZ 

z 


.1406 
"  T764- 

=  .0857  m3 
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SAMPLE  CALCULATION  #3 
The  following  load  calculations  are  examples  using 
values  taken  from  Run  III  of  the  hogging  case. 

NOMENCLATURE  AND  VALUES 


T  =  draft 

=  15»2  meters 

P    =   density  of  water 
=   1.026  mtons/nr 

M  =  moment  at  frame  195 
=  484,900  mtons-meter 

Q  =  shear  at  frame  195 
=  610  mtons 

W  =  static  steel  weight  of  quarter  tank  section 
=  74.55  tons/meter 

a  =  acceleration  of  gravity 
S         /   2 
=  9.81  m/sec 

a,  =  acceleration  due  to  heave 
h  2 

=  2.76  m/sec 

a  =  acceleration  due  to  pitching 
P  2 

=  . 04  m/sec 

n  =  number  of  joints  in  beam  model 
=  208  joints 

c   =  static  load  of  cargo  (quarter  tank  section) 
=  3607  mtons 

r.  =  ratio  of  actual  ship's  bean  to  model's  beam 
=  (22,00)  =  1  Q5 
v20.78;    1<u-» 

r-.  =  ratio  of  actual  quarter  tank  section  length  to  model's 
1        length 

=  (21.7/21.06)  =  1.030 
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Forces  Due  to  Water  Acting  on  the  Hull 


P  =  pressure  on  bottom  hull  =  (T)(/o) 
=  (1.026)(15.2) 

=   15*6  intons/cm 

The  water  pressure  loads  were  applied  to  the  vertical 
members  as  indicated  below. 


ISA  -m 


m  =  (P) (member  width) 

=  (15.6M2.55) 
=  39.78  mtons/cm 

F  for  the  vertical  members  at  frame  228  is  half  of  the 
m 

above  value  due  to  symmetry  conditions. 


The  loads  on  the  bottom  were  divided  such  that  half  of  the 
loads  were  applied  to  longitudinal  members  and  half  were 
applied  to  the  transverse  members. 

F  =  (P) (member  width) (r,  ) 

=  (15.6)(2. 55)11.1059) 
=  21.06  mtons/m 

F  for  transverse  members  1,  10,  19,  28,  38,  49,  and  60  is 

W 

one  half  the  above  value  due  to  symmetry. 

The  forces  on  the  longitudinal  member  were  calculated  in  a 
similar  fashion  using 


F-l  =  (P)  (member  width)  (r^ 
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Loads  Due  To  Cargo 

g  =  multiple  of  force  felt  by  the  quarter  tank  due  to 
dynamic  conditions 

s  ( g    -j-  3    4'  ?!   )  /  » 

o       '• l       iJ       o 

=  (9,81  1+  2.?6  -t"  .04)/  9,81 
=  1.4 


F  =  dynamic  cargo  load  ~  (c)(g) 
=  (3607) (104) 
=  50.50  mtons 


There  are  11  support  joints,  two  of  which  only  carry  half  of 
the  load  of  the  others,  due  to  symmetry.   Thus* 

f  =  5050/(9  +  ,  5  J-  ,  5) 

=  50,5  mtons/support  joint 

Loads  Due  To  Steel  Weight 

F  =  ( w )  (  a ) 

w    ■  '■'-'■ 

=  (74,55)(21  .7)0  <'0 

=  2?.65  mtons 

f  ~  F-,  /  n 

w    w 

=  2265/  208 

-  10,89  mtons/  joint 
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BOUNDARY  LOADS  -  FRAME  195 

The  horizontal  neutral  axis  was  calculated  to  be  located 
10.5  meters  above  the  bottom  deck.   The  loads  simulating 
the  bending  mement  were  assumed  to  be  distributed  as  indi- 
cated in  Figure  30.   The  moment  that  had  to  be  applied  to 
the  quarter  tank  model  was  one  half  of  the  actual  moment  of 
the  whole  ship  at  fr-ame  195«   The  forces  applied  above  the 
neutral  axis  would  form  half  of  the  couple.   Thus  to  deter- 
mine F.  and  F, 


(2)(2)  =  1/2  couPle 

F.   Lp    ? 

♦  (f)cfH-T"L2) 

where      L1  =  distance  between  the  inner  and  outer  webs 
=  8.27  m 

L  =  distance  between  the  neutral  axis  and  the 
main  deck 

=  13. ^5  m 

F  F  F 

U8U^°     =  ^(8.37)(13.4)  +  -^(8.37)(13.*5)  +  ~J   (13.^5)2+ 

t  yd3.^5)2 


F,  ~  779.9^  mtons/meter 
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1=   Fb(L3)(L4)  +  (2)^  (|L^) 

where      L„  =  width  of  double  bottom 

=  20.78  m 

L^  =  distance  from  double  bottom  to  neutral  axis 
=  10.5 

^8y°   =  Fb  (20.78)(10.5)  +  Fb  (10. 5)2  (2/3) 
F,  =  415.6  mtons/meter 

The  shear  load  at  frame  195  for  the  quarter  tank  section 
was  610.5  mtons.   This  was  distributed  equally  between  the 
ten  joints  shown  in  Figure  30*   Thus 


q  =  Q/1.0 

=   6l1°0'  *?  =   61.05  mtons 


Boundary  Loads  -  Frame  228 

In  order  to  calculate  the  shear  forces  at  frame  228  to  pre- 
vent rigid  body  motion  of  the  model  it  is  necessary  to  make 
the  summation  of  the  forces  in  the  vertical  direction  equal 
to  zero. 

F  =  summation  of  Z  direction  forces 

z 

=  Q  -  F   -  F  +  water  on  bottom  hull 

=  610.5  -  2265  -  5050  +  7^37 
=  742  mtons 
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F 
[228  10 

242 

10 


-  7^.2  mtons/ joint 
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APPENDIX  B 
Computer  Input  for  Run  Iir 
of  Hogging  Case 
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